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Abstract: A nonperturbative procedure of solving the time-dependent 
Schrodinger equation, called the multi-projection approach or phase dy- 
namics of quantum mechanics, is derived and illustrated. In addition to 
introducing a method with that time-dependent systems become solvable 
(under the assumption that corresponding time-independent systems are 
solvable), the new approach unveils several misconcepts related to the usual 
wavefunction expansion and the standard perturbation theory. 

PACS numbers: 03.65-w 



1 Introduction 

The standard methodology of quantum dynamics, initially due to Dirac[|]], 
assumes that a wave function can always be expanded into a series of eigen- 
functions 

^C,(t)e-^^(r), (1) 

where Wfc(r) represent a set of eigenf unctions, which is, according to the 
customary understanding, associated with the initial Hamiltonian H(to) of 
the quantum system. By inserting the expansion into the time-dependent 
Schrodinger equation and doing some mathematical manipulations, a set of 
coupled differential equations is obtained as 

ih dCk = j2c m e i ^-^ t V km , (2) 

at m 

where Vtm are the matrix elements of the Hamiltonian variation V(t) = 
H(t) — H(to). A common concept, taught everywhere and established solidly 
in our mind, is that the expansion (|l]) is the formal solution of the wave func- 
tion and the equation set (|2|) describes the system's behavior as exactly as 
the original Schrodinger equation does. Many perturbative and nonpertur- 
bative approaches are based on or related to this concept. (As an example, 
see the derivation of the Liouville theorem in quantum statistics.) 

Our studies, however, showed that the concept outlined above, while 
enjoying tremendous successes, involved serious difficulties [Q]|3|. Some of 
them are the following. 

• In numerical computations, the sum of squared moduli of all coeffi- 
cients, i.e., J2\Ck(t)\ 2 , can easily go to infinity if V(t) is not truly 
small. That is to say, the Dirac theory, like many other types of per- 
turbation theories, suffers from the divergence difficulty. 

• The notion associated with (|l|) and (||) asserts that if a quantum sys- 
tem, say a harmonic oscillator, is initially in the ground state, it will 
make transition, partially though, to a higher state at the next mo- 
ment, and then to a higher higher state at the next next moment. The 
dynamical process will never terminate if H (t) — H(to) is nonzero as 
t — > oo (another symptom of the divergence). In contrast with that, 
a general analysis of quantum mechanics states that whenever the 
Hamiltonian H(t) becomes steady, equal to H(to) or not, the system 
will instantly be settled in a stationary state. 



• According to the state-transition picture provided by the Dirac the- 
ory, a ground-state system will never lose its energy under any circum- 
stances. This conclusion does not appear to be in harmony with other 
theoretical and experimental observations. 

• The fact that an electromagnetic field can be represented by different 
gauge fields implies that there are an infinite number of equation sets 
taking the form (|2|) for one definite disturbance. Direct numerical 
calculations show that these equation sets are not equivalent to each 
other. 

Many questions then arise, of which some must be of fundamental inter- 
est to the community. How does the Dirac formalism, derived mathemati- 
cally from the Schrodinger equation, involve so many difficulties? Why can 
the theory, running into problems ultimately, still offer good results under 
certain conditions? What on earth are these conditions? Is there a method 
by that the time-dependent Schrodinger equation can be solved more ad- 
equately? Among them, the last question seems to be the most essential 
one. If we find out a proper way of solving this fundamental equation, many 
questions, in particular those related to the Dirac perturbation theory, can 
be answered accordingly. 

Although a different perturbation theory[||, which corresponded to a 
similar theory in classical mechanics Hj and suffered less problems, was put 
forward, the issue was far from cleared-up and the aforementioned questions 
remained to be open. 

In this paper we wish to report on a nonperturbative procedure of solving 
the time-dependent Schrodinger equation. Interestingly, the new procedure, 
called the multi-projection approach or phase dynamics of quantum mechan- 
ics, is still based on Dirac's idea: solving a nonstationary quantum system 
with knowledge of stationary quantum systems. But, unlike its predecessor, 
this approach suffers from no divergence difficulty and exhibits surprising 
effectiveness. Any dynamical system, subject to weak or strong fields, be- 
comes analyzable and calculable provided that the corresponding stationary 
states can be solved analytically or numerically. 

The structure of this paper is the following. Section 2 introduces our new 
approach, which is based on an assumption that the time-dependent Hamil- 
tonian can be approximated by its stepwise varying counterparts. In Sec. 
3, similarities and dissimilarities of the proposed approach to the influen- 
tial path-integral approach are remarked. Simple applications are presented 



in Section 4, where effectiveness of the new approach is illustrated. Sec- 
tion 5 rederives the Dirac perturbation theory and answers several related 
questions. Sec. 6 gives discussions on the common concept of wavefunction 
expansion. Sec. 7 concludes the paper. 

2 Multi-projection approach 

In this section, we try to determine the wave function ^(t) on the condition 
that the initial wave function \P(to) and the time-dependent Hamiltonian 
H(t) are given. 

According to the basic formalism of quantum mechanics, the solution of 
the Schrodinger equation can be written as 

rt 

y(t) = e~ iJt o H{r)dT y(t ). (3) 

The solution above is just a formal one except in the situation where the 
involved Hamiltonian is independent of time. To generally evaluate it, we 
slice the time interval from to to t into n segments, equal to each other or 
not, as 

Ati=ti-t , •••, •••, At n = t-t n - 1 . (4) 

The formal solution (0) can then be rewritten as 

-if* HMdr -1 f' 1 HMdr , „ 

Or, in terms of the intermediate quantum states, 

*(t J -) = e~* J «/-i H(T) *'*(t i _ 1 ) (J = 1,2,---). (6) 

Referring to Fig. 1, we explore the possibility to replace the Hamiltonian 
H(t) with its stepwise varying approximation H{t). Without losing gener- 
ality, consider a charged particle in a time-dependent electromagnetic field, 
whose Hamiltonian reads 

tf(*) = ^(p-^A) 2 + <I>. (7) 

2m c 

It is rather obvious that if <3? ^ and A = 0, the replacement of H{t) 
with H(t) is justified by the observation that the two Hamiltonians rep- 
resent roughly the same physical system. For the case in that both the 



scalar potential 3> and the vector potential A are nonzero, the situation 
becomes somewhat complicated. It is well-known, the vector potential A 
can be separated into two parts: the longitudinal field and the transverse 
field p]. By realizing that the longitudinal vector field can cause trouble to 
our formalism (the reason for saying it will be clear), only the scalar field 
and transverse vector field will be treated herein. This should be allowable 
in view of that an appropriate gauge transformation can always make the 
longitudinal vector field vanish. Under these understandings, we define the 
Hamiltonian H(t) as 

Hj(t) = 4r f° H(t)dt for {tj-x < t < tj); (8) 

or, even in a simpler way, Hj(t) = [H(tj-x) + H(tj)]/2. Thus, during each 
of the time segments defined by (H) the new Hamiltonian Hj is independent 
of time and the typical intermediate state in @ becomes 

Jf^we-i^^fe-i). (9) 

The problem is reduced to a familiar one related to solving stationary sys- 
tems. For each of the intermediate Hilbert space associated with Hj, we 
have normalized intermediate energy eigenfunctions 

H J wi(r)=Elwi(r). (10) 

After the wave function ^f(tj-x) is known, the wave function \P(tj) can be 
expressed by 

*(*;) = £^e-K At ^(r), (11) 

where u k = E J k /h and C° k is determined by a projection 

C{ = f ntj-x)Wl\v)dv. (12) 

Or, in the Dirac notation, 

\t J )=J2(k J \t j -x)e-^ At >\k J ), (13) 



where \kj) stands for Wl. Eq. (13) shows that for a short time interval 



a dynamical system and its correponding stationary system evolve in the 



same way. By repeating the step as presented above, the wave function at 
the time t can be expressed as 






I*) = >:•■• >J e- l ^^ Atj \k n )(k n \k n -i) ■ ■ ■ (h\t ). (14) 



K>n tZ\ 



For convenience of discussion, we will name \k n ) (k n \k n -i) ■ ■ ■ {ki\to) as a 
minute component of \to), by which we refer to the fact that such a compo- 
nent is in usual situations very small. It is easy to prove that the sum of all 
the minute components, with no phase factors involved, is just equal to the 
initial wavefunction. 

In deriving the above formalism, it has been assumed that the wave 
function ty(t < to) does not contain any phase factor in the form e l ^ t ' r ^ 
where / depends on t and r explicitly. Generally speaking, e l ^ t,r > is not a 
uniformly continuous function with respect to t and r and great difficulties 
will be encountered if we try to expand a wavefunction having such phase 
factor. A related and important point is that since we have assumed that all 
Hj contain no longitudinal vector fields, all the intermediate wave functions 
are also free from such phase factors. 

Numerically speaking, the accuracy of this method largely depends on 
the step size. If all steps in a calculation approach infinitesimal ones, the 
resultant wavefunction will approach the real dynamical one. 

For a dynamical system, the energy at any moment can naturally be 
defined by the intermediate Hamiltonian and the intermediate wavefunction. 
In the sense of taking limit, this definition is an accurate one. 

Before finishing this section, we wish to point out that this approach of- 
fers a different picture about quantum dynamics. The Dirac theory leads one 
to imagine how a nonstationary system makes transition from one eigenstate 
to another. In contrast with that, the multi-projection approach reveals 
that a nonstationary system is subjected to no other change than that each 
minute component of the wavefunction acquires its own dynamical phase 
factor. To be consistent with the spirit of enlightening discussions about 
Berry's phase||, the obtained formalism may be called phase dynamics of 
quantum mechanics. 

3 Equivalence to the path-integral approach 

In this section, the path-integral approach [0]|§| is briefly reviewed and then 
the relation between the path-integral approach and the multi-projection 
approach is remarked. 



In the path-integral approach we have 

*(x,t) = (x,t\xo,to)y(xo,to)dxo, (15) 

where (x, t\xo, to) is named as the propagator. (For simplicity, only the one- 
dimensional case is considered.) It is found that the propagator takes the 
form 

(x,t\xo,to)= f e lS[x{T)]/n D[x(r)i (16) 

where Z)[x(r)] is the measure associated with all possible paths x{t) and 
the action S is the integral of the Lagrangian along the path 



S 



in 



L(r)dT 



t r 



u> 



-mv — V(x(t),t) 



dr. 



(17) 



The formalism above is quite formal unless the discrete form of it is written 
out. By slicing the time (t — to) into N equal segments and denoting (t — 
to)/N by e, the action integral becomes 



JV-l 
j=0 L 



m{x j+ i -Xj) 
2e 



eV(x',e' 



(18) 



where x' and e' take values within the intervals (xj,xj+i) and (tj,tj+x) 
respectively. Upon this, the propogator becomes 



^oo +oo 



lim A 

N^oo 



-oo — oo 



JV-l ■ 

ex P £ I 

j=0 



m{Xj + i — Xjj 

2e 



eV(x',e') 



dx\--dxN-\i (19) 



where the factor A can be determined by the limit (x,e|xo,0) — » 5{x — Xq). 
For purposes of this paper, we wish to take a look at how the path- 
integral formalism yields the standard Schrodinger equation ||. Consider 
the propagator associated with one slice of time after t = 



/ m \V 2 ( i 

which means 



m(x — xqY 

2e 



eV(x',e' 



(20) 



t, ^ ( m Y /2 



, exp u 



m(x — xo) z 

2e 



eV(a/,e?) 



V(x o ,0)dx . (21) 



By noting that the region (x — xq) 2 ~ 2ehir/m gives the main contribution 
to the integral, we obtain, after taking approximations, 



$)(x,e) - *(x,0) 



2m dx 2 



9(x,Q), (22) 



which is just the finite differential form of the Schrodinger equation. 

The similarities between the path-integral and the multi-projection are 
found by the following observations, (i) Both approaches acquire clear mean- 
ing in their discrete forms, (ii) Both approaches are approximate theories 
in their discrete forms and convergent to the exact theory as steps of the 
discrete forms become infinitesimal, (iii) For a short time interval, both 
approaches show indifference to whether or not the quantum system is truly 
time-dependent. In the path-integral approach outlined in this section, if 
we replace the Lagragian with its stepwise varying approximation, all the 
formulas, in particular (p0|), (p|) and (p2]), make no essential change; as the 
steps go smaller and smaller, exactly the same limits will be obtained. All 
these imply that the equivalence of these two approaches is indeed there. 

It is now worth mentioning dissimilarities. The path-integral approach 
formulates the time unitary transformation of wave function in the ordinary 
spatial space, whereas the multi-projection approach does the same job be- 
tween a series of Hilbert spaces whose bases are formed by eigenfunctions of 
the intermediate Hamiltonians. Due to this difference, the formalism in the 
path-integral approach takes a form of probability integral, while the for- 
malism in the multi-projection approach takes a form of a series of definite 
projections. In terms of studying practical systems, it is expected that the 
two provide different computational ease. 

Finally, a note about gauge. It is rather necessary for both the ap- 
proaches to assume no longitudinal vector field to exist. 

4 Simple applications 

Some concrete examples, in which effectiveness of the multi-projection ap- 
proach manifests itself, are presented in this section. 

Firstly, consider a case in that the potential of a harmonic oscillator is 
subject to a "sudden" change and the Hamiltonian reads 

H(t) = ^ + ^kx 2 , (23) 



with 



*«>=;!!; 8, 



in which r) is a positive constant other than 1. Physical intuition says that, 
if r) > 1 the total energy of the system must go larger (the spatial room for 
the oscillator is "compressed"); if r] < 1 the total energy must go smaller. 

Though the situation given above is terribly simple, it still defines a dy- 
namical system that needs to be solved in one way or another. As indicated 
in the introduction, if the standard method is applied serious difficulties 
arise. In what follows we try to use our multi-projection approach. 

Assume that the system is initially in the ground state, namely|| 

*(0) = N (a) exp(-aV/2), (25) 

and the wave function after t = takes the form 

*(*) = J2C n e- iin+ ^' t N n (a')H n (a'x)e- 1 2 a ' 2x2 . (26) 



The notation in Eqs. ( g5J ) and (26) is rather standard: a = (mk/H ) 1 ' 4 , 
a' = tfrja, uj 1 = ^/rjuJ = \frjkjm and 

/ a \ 2 
N n (a) = (j^lj , H n+1 (0 = 2ZH n - H n ^ (with H = 1). 

The energy of the system at t = is known to be E$ = 0.5hco. By making 
use of (]T2|), we obtain, for r\ = 0.5 2 , 

i i 

4 2 / 1 \l 



Co=[- g ) , C 1 = 0, C 2 = --{-) , C 3 = 0.... (27) 

Numerically, we have 

C « 0.9710, C 2 « -0.2289, C 4 « 0.0661, C 6 « -0.0201; (28) 

these imply ^ |C n | 2 ~ 1 and the final energy becomes 

(E) = J2 \Cn\ 2 (n + 0.5) W « 0.6246^0- (29) 

Similarly, we obtain, for rj = 0.9 2 and rj = l.l 2 respectively, 

(E) ss 0.9050£ and (£) « 1.1050£ - (30) 



It is now obvious that with the multi-projection approach (i) the normal- 
ization condition for wave function, namely J2 \C n \ 2 = 1, is automatically 
satisfied; (ii) whenever the Hamiltonian becomes steady the system is set- 
tled in a stationary state; (iii) even if a quantum system is initially in a 
ground state, the total energy of it can decrease or increase; (iv) no gauge 
issues pose problems if we let the longitudinal vector potential vanish al- 
ways. In other words, all the difficulties troubling the standard dynamical 
theory disappear completely. 

The following similar example shows that the proposed method is a good 
mean for evaluating the phase shift of a quantum system. If the oscillator 
expressed by ( ]23|) is again in the ground state at t = and the function S(t) 
in (p5|) is, instead of (|24|), 

r i (t < o) 

S(t) = I r) (0 < t < T) (31) 

[ 1 (t > T), 

the formulation of the above example can still be employed except that a 
projection at t = T onto the original Hilbert space is needed. In general, the 
final wave function involves energy transition and phase shift. As a special 
case, we consider the situation in which T = 4tt/uj' . Simple calculation tells 
us that at t = T the wave function will come back to its original value and 
original phase (by virtue of e l27T = 1) as if the system has been completely 
"frozen". But, on the other hand, one finds that if the system gets no 
disturbed, it will during the same time acquire the phase factor 

exp(-iwr/2) = exp(-i27r/v^). (32) 



It is then obvious that the disturbance (|3l]) makes the wavefunction have 
the additional phase factor exp(i2n / ^/rj) . 

An inspection of the last example states that if a perturbation is capa- 
ble of changing the system's eigenfrequencies, it can in general make the 
wavefunction have an additional phase shift, irrespective of whether or not 
there is energy transition. As shown by experiments flO|| , such phase shift 
has physical effects and should be treated with care. 

5 Rederivation of the Dirac perturbation theory 

In this section we rederive the Dirac perturbation theory with help of our 
proposed approach. The purpose of doing that is two folds. One is to 
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illustrate the analytical ability of the approach; and the other is to answer 
questions related to the Dirac perturbation theory. 

To obtain an analytical formalism comparable to the Dirac one, we pre- 
suppose that (i) the quantum system of interest is initially in the nth eigen- 
state of Hq, the perturbation applies at t = and vanishes completely at 
t = T; (ii) the perturbation can be approximated by a series of pulses whose 
values rise sharply and vanish sharply, as illustrated in Fig. 2 (the pulse-like 
perturbation has the same physical effects as the real perturbation); (iii) the 
system's wavefunction involves no additional phase-factor shift before and 
after each of the pulses. With these assumptions adopted, the derivation 
here can be deemed as a simple application of the "sudden approximation" 
due to Pauli(ll|] 



Consider one specific pulse of the perturbation that exists between t' and 
t' + At'. The leading term of the wave function at t = t' is e~ tuJnt |n), at 
t = T it becomes, by repeated use of dl3|) , 

e- iuJnt '\n) -^J2 e ~ iUlnt ' e ~ iUlkAt '( k \ n )\ k ) 

k 
-^J2 e ~ iUJ " t ' e ~ iUJkAt ' e ~ iUlm{T ^'' At ' ) (m\k)(k\n)\m), (33) 

k,m 

where \m) stand for eigenfunctions defined by the Hamiltonian Hq and \k) 
eigenfunctions defined by the intermediate Hamiltonian within the pulse. 
Note that in ( p3[) all the phase factors take on their unperturbed values 
before t' and after t' + At' as if no other pulses exist, which should in general 
be regarded as a rough approximation, as revealed in the last section. Since 
At' is short, the following approximation is acceptable 

]T(m|£;)e-^-^ At ' (k\n) « (m 1 - ^W(t') - H ]At' n) , (34) 

k \ n / 

thus, the rath coefficient of the wave function at t = T is 

b m = —At' {m\V\n)e- i{ -^-° Jm)t ' &- itj3mT , (m ± n) (35) 

where V = H(t') — Hq. Taking contributions from all pulses into account, 
we obtain at t = T 

b m = - % -e- l ^ T [ T Vrnne-^-^'dt' (m^n). (36) 

n Jo 
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Except the phase factor exp(— iu> m T), the formula above is consistent with 
the well-known one. 

Unlike the standard derivation in textbooks, the present derivation clears 
up several subtle and important things, of which some were addressed in the 
introduction. Firstly, it indicates that if the final Hamiltonian H(t > T) is 
not the same as the initial Hamiltonian, this formalism will assume that 
"perturbation pulses" are always there and thus yield misleading results. 
Secondly, it stresses that if the perturbation, represented by the Hamiltonian 
variation V(t), is not relatively small or if the action time of the perturbation 
is not relatively short, the leading term of the wavefunction, including its 
phase factor, will in general be disturbed significantly and the accuracy of 
the formalism will partly or entirely be demolished. Thirdly, it suggests that 
higher-order solutions directly given by the Dirac perturbation theory are 
not meaningful. Finally, it implies that in order for the formalism to hold, 
an appropriate gauge, like the one defined in Sec. 2, needs to be adopted. 

6 Discussion on the wavefunction expansion 

Let us return to the subject put forth at the very beginning. What have 
been presented in this paper, together with those in Ref. 2, actually bring 
out that the usual concept about the wavefunction expansion bears problem- 
atic aspects. The formulation in Sec. 2 clearly shows that after a quantum 
system leaves its initial state, the eigenfunctions and eigenfrequencies asso- 
ciated with the initial Hamiltonian become out-of-date. If we forcefully use 
them to express the wavefunction at the later times, the coefficients of the 
would-be series have to adjust themselves violently — so violently that they 
cannot be regarded as continuous variables in the usual sense. If the aim 
of a calculation is to determine these coefficients continuously, the situation 
will soon become out of control: the more terms are taken into account, the 
faster temporal scales and the larger spatial scales get involved, the bigger 
and more serious errors thus enter the scheme. 

In other words, although one can, at any fixed moment, expand a dy- 
namical wavefunction into a series of the initial eigenfunctions as, up to a 
phase factor, 

J2C n (t)W n (r) or Y, C ^ t )e~ iuJnt W n (r), (37) 

it is not appropriate to assume that a dynamical wavefunction can be repre- 



seted by an expansion that takes the form (37) and has continuously time- 
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varying coefficients. 

The arguments that have just been presented unveil, in another perspec- 
tive, that the concepts introduced by this approach — intermediate Hilbert 
spaces, intermediate eigenfrequencies and discrete projections between in- 
termediate Hilbert spaces — are more fundamental than they appear to be. 
(Also, it justifies the path-integral approach, whose objective is to formulate 
a wavefunction at certain discrete moments.) 

7 Summary 

We have proposed a nonperturbative procedure, called the multi-projection 
approach, which in a way unifies treatments of different quantum systems: 
stationary and nonstationary, weakly-disturbed and strongly-disturbed. In 
view of that knowledge about stationary systems has been accumulated for 
long and many systems with strong fields need to be studied, such unification 
appears to be quite desirable. 

The equivalence of the multi-projection approach to the path-integral 
approach also suggests that the proposed method may find its applications 
in a variety of quantum fields, though it is still too early to say what exactly 
they are. 

This paper has shown that a dynamical process of quantum system is, 
rather strikingly, characterized by the phase-factor evolution of each minute 
component of the wavefunction. This "state transition" picture is quite 
different from the standard one provided by the Dirac perturbation theory. 
The standard picture, associated usually with one fixed Hilbert space and 
paying much less attention to disturbance of phase factors, is applicable only 
if some conditions are simultaneously satisfied, which in general include: the 
final Hamiltonian of the system is the same as the initial Hamiltonian, the 
perturbation is relatively small, and the action time of the perturbation is 
relatively short. 

We have pointed out that a dynamical wavefunction cannot be repre- 
sented by a wavefunction expansion that has continuously time-varying co- 
efficients. 

Stimulating discussion with Professor Han-ying Guo is gratefully ac- 
knowledged. This paper is partly supported by School of Science, BUAA, 
PRC. 
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Fig. 1, A typical time-dependent Hamiltonian and its stepwise- varying ap- 
proximation. 



,V(t) = H(t) - H 



1 



Fig. 2, A Hamiltonian variation and its pulse-like approximation. 
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